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Abstract 

In this paper we study lattice rules which are cubature formulae to approximate 
integrands over the unit cube [0, l] s from a weighted reproducing kernel Hilbert 
space. We assume that the weights are independent random variables with a given 
mean and variance for two reasons stemming from practical applications: (i) It is 
usually not known in practice how to choose the weights. Thus by assuming that 
the weights are random variables, we obtain robust constructions (with respect to 
the weights) of lattice rules. This, to some extend, removes the necessity to carefully 
choose the weights, (ii) In practice it is convenient to use the same lattice rule for 
many different integrands. The best choice of weights for each integrand may vary 
to some degree, hence considering the weights random variables does justice to how 
lattice rules are used in applications. 

In this paper the worst-case error is therefore a random variable depending on 
random weights. We show how one can construct lattice rules which perform well for 
weights taken from a set with large measure. Such lattice rules are therefore robust 
with respect to certain changes in the weights. The construction algorithm uses the 
component-by-component (cbc) idea based on two criteria, one using the mean of 
the worst case error and the second criterion using a bound on the variance of the 
worst-case error. We call the new algorithm the cbc2c (component-by-component 
with 2 constraints) algorithm. 

We also study a generalized version which uses r constraints which we call the 
cbcrc (component-by-component with r constraints) algorithm. We show that lat- 
tice rules generated by the cbcrc algorithm simultaneously work well for all weights 
in a subspace spanned by the chosen weights 7^, . . . , 7 ( - r * ) . Thus, in applications, 
instead of finding one set of weights, it is enough to find an r dimensional convex 
polytope in which the optimal weights lie. The price for this method is a factor r in 
the upper bound on the error and in the construction cost of the lattice rule. Thus 
the burden of determining one set of weights very precisely can be shifted to the 
construction of good lattice rules. 

Numerical results indicate the benefit of using the cbc2c algorithm for certain 
choices of weights. 
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1 Introduction 



In this paper we study the integration error using a lattice rule. A lattice rules is a 
quadrature rule of the form 

<w/>4x:'({m 

n=0 

where N > 1 and s are natural numbers, g G {1,...,N — 1} S , and where {x} = 
({xi}, . . . , for a; = (x\, . . . ,x s ) and {x} = x — [_xj stands for the fractional part of 

a nonnegative real number x. Lattice rules are quasi-Monte Carlo algorithms which are 
useful to approximate integrals Jj ^ f(x) dx. In the following we present a survey of the 
literature and some background on lattice rules. 

1.1 Literature survey and background on lattice rules 

It has been shown that lattice rules are efficient for approximating integrals of periodic 
functions (see [HI [Hll [El El])- The construction of lattice rules has seen many advances in 
recent years. An important framework in which to study lattice rules (and other quadra- 
ture rules) are reproducing kernel Hilbert spaces, which have first been considered in 
[7] and are now a standard tool in quasi-Monte Carlo integration. The component-by- 
component (cbc) construction (where the generating vector g is constructed one compo- 
nent at a time) was first discovered by Korobov [10J and P, Theorem 18, p. 120] and 
independently rediscovered by Sloan and Reztsov |24j . In [221 122] this idea has been fur- 
ther developed to allow one to use lattice rules also for nonperiodic integrands. Optimal 
convergence rates for lattice rules constructed this way have been shown in [HI HO] and 
independently in [TTJ for a prime number of points iV and in [5] for a nonprime number of 
points. A breakthrough in reducing the construction cost of the cbc construction has been 
achieved by Nuyens and Cools in [161 Ej> who showed how the fast Fourier transform 
can be used to reduce the construction cost of the search algorithm. A further very im- 
portant development has been the introduction of weighted function spaces by Sloan and 
Wozniakowski [25] . Therein, the authors make the important observation that integrands 
may have different dependence on different projections. To take this fact into account, the 
authors introduced so-called weighted function spaces which yields a weighted worst-case 
error criterion. A comprehensive introduction to weighted function spaces and tractability 
questions as well as further background can be found in the comprehensive monographs 

[IHII51. 

1.2 Worst-case error 

In the following we introduce a specific reproducing kernel which will be sufficient to 
illustrate our algorithm. In order to keep the notation as simple as possible, we do not 
consider the most general case possible. 

Consider the reproducing kernel (see pQ) K : [0, l] 2 — > C defined by 

Kj(x,y) = ^B 2 {{x-y}), 
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where 7 > is a nonnegative real number, the 'weight', and ^(z) 



z 



z + 1 /6 is the 



Bernoulli polynomial of degree two (cf. [7]). The Bernoulli polynomial B 2 has the Fourier 
series 

fcez\{o} 

The reproducing kernel K 1 defines a reproducing kernel Hilbert space % 1 of absolutely 
continuous, periodic functions on [0, 1] which integrate to 0, with inner product 



(f,g) 



2^ 

7 



E fc2 /(^)- 



kez\{o} 



For dimensions s > 1 we consider the reproducing kernel 



KJx,y) 



E ^Y[B 2 ({xi-y t }), 



where <S = {1, . . . , s}, 7 = (■y u )fD^ucs is a set of nonnegative real numbers 7 U associated 
with the projection onto the coordinates in u (we refer to these numbers as the 'weights' 
[25]), ) and y = (y±, . . . ,y s ). The associated reproducing kernel Hilbert 

space is denoted by "H 7 which is a sum of tensor products of the reproducing kernel 
Hilbert space with kernel K 1 and the space of constant functions, see [6j [7J |25] for more 
information. The kernel K n can also be interpreted as the shift-invariant kernel of a 
reproducing kernel Hilbert space of non-periodic functions [SJ. Thus the results here 
can also be interpreted for randomly shifted lattice rules in the associated non-periodic 
reproducing kernel Hilbert space (as for instance in [XT], IT6] ITT] and many other papers). 

The integration error using a lattice rule with generating vector g e {1, . . . , N — 1} S 
is defined as 



e{U- v P N (g))= sup 

feH~,,P N (g) 
II/II«,<1 



[0,1]' 



n=0 



It was shown in (6j [7] that 

e 2 (H- r ,P N (g)) = / / K 1 (x,y)dxdy 

J[o,i] s J[o,i] s 

N-l 



N-l . N-l 

~ N E / A ^' ^({ng/N}, {n'g 

n=0 ^[O, 1 ]" n,n'=0 



/N}) 



N-l 



e 4^nw}), 



where the last inequality follows from L K 7 (x, y) dx = 1 and the fact that (see P UJ) 

1 JV-1 1 JV-1 

- E K 1 {{ng/N},{n'g/N}) = -Y,K 1 {{ng/N},Q). 



N 2 



n,n'=0 



n=0 
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We use the last expression as error criterion in this paper, i.e. 

e 2 (^ 7 ,/M<7))= £ (lu^llBittnm/N})). (1) 

As indicated above, ([I]) can be interpreted as the square worst-case error in a Korobov 
space or the mean square worst-case error of a randomly shifted lattice rules in a Sobolev 
space. The square worst-case error e 2 ("H, Pjv(gr)) is commonly used as error-criterion in a 
cbc construction EH UHl H3 E21 E3] . 

The cbc algorithm is a greedy search algorithm to find a good generating vector 
9* = idi j • • • 5 5s ) £ { 1 j • • • ) N — 1 } . This algorithm works the following way. First one 
chooses a number of points N, the dimension s and some weights 7. The component-by- 
component construction then finds a generating vector g* = [g* , . . . , g*) in the following 
way: 

• Set g{ = 1. 

• For j = 1, . . . , s set 

g* = fagmm. 1 < e < N _ 1 e 2 (?lj, P N ((gt, <£_ x , 2))). 
1.3 The aim of the paper 

In practice, it is usually not known how to choose the weights j u , 7^ u C 5 in the 
worst-case error criterion. Some suggestions on how to choose the weights in financial 
applications have been put forward for instance in [26l [27] . Another method of choosing 
the weights is by choosing them such that the error bound is minimized [T2]. However, 
choosing good weights remains a particular challenge for the application of lattice rules. 

In this paper we assume that the weights are independent random variables with a 
given mean and variance for two reasons stemming from practical applications: 

• It is usually not known in practice how to choose the weights precisely. By assuming 
randomness in the weights permits a 'measurement error' or noise in choosing the 
weights. 

• It is convenient to use the same lattice rule for many different integrands. The best 
choice of weights for each integrand may vary to some degree, hence considering the 
weights random variables seems to be the right model in this case. 

Indeed, it is desirable to have quadrature rules (lattice rules) which are robust with 
respect to the weights (7^), that is, for which one obtains a good convergence behavior, 
not only for one given choice of weights, but for a whole range of weights. In order 
to construct lattice rules which have this property, we assume that the weights are not 
given (or fixed), but rather, we assume they are chosen randomly with a given mean and 
variance. In this way, the square worst-case error e 2 ("H 7 , P^ig)) is a random variable 
(with respect to the weights 7). In the following we propose an algorithm to construct 
lattice rules which have a small expectation value and, at the same time, a small variance 
of e 2 {7i-y, Pjv(fl')) with respect to the random choices 7. In a nutshell, the existence 
of such a lattice rule is guaranteed by the fact that more than half of the generating 
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vectors have small expectation value of the worst-case error and more than half of the 
generating vectors have small variance of the worst-case error. Thus there exists at least 
one vector for which both, the expectation value and the variance, are small. This yields a 
component-by-component algorithm with 2 constrains, which we call the cbc2c algorithm. 

We also study a general version which uses r constraints which we call the cbcrc 
(component-by-component (with) r constraints) algorithm. We show that lattice rules 
generated by the cbcrc algorithm simultaneously work well for all weights in a subspace 
spanned by the chosen weights 7^, . . . , j^. Thus, in applications, instead of finding one 
set of weights, it is enough to find an r dimensional convex polytope in which the optimal 
weights lie. The price for this method is a factor r in the upper bound on the error and 
the construction cost of the lattice rule. Thus the burden of finding one set of weights 
can be shifted to the construction of good lattice rules. 

Theoretically one could make an exhaustive search to obtain a lattice rule which 
simultaneously works well for all choices of weights. This may eventually shift the question 
of how to choose the weights for a particular problem to the problem of finding a universal 
lattice rule which simultaneously works well for all choices of weights, thereby removing 
the need to choose weights in the first place. The computational challenge though is, that 
for higher dimensions finding such a lattice rule is currently intractable (since the cost 
depends exponentially on the dimension). Further, also the upper bound from this paper 
depends exponentially on the dimension when r = 2 s — 1. This method may be useful 
though for integrands with low truncation dimension d by choosing r = 2 d — 1 in this 
case. 

We note that a similar theory can be applied to polynomial lattice rules and related 
point sets [3J US] • 

In the next section we study the implications of the assumption that the weights are 
random on the square worst-case error. In Section [3] we first repeat some important 
insights from [161 E]- We introduce the cbc2 algorithm and show that the constructed 
lattice rules work well for weights taken from a set of large measure. We then consider 
the cbcr algorithm and consider the geometrical interpretation of the algorithm. It is 
shown that the square worst case error satisfies a certain bound for all weights in an r- 
dimensional convex polytope which is defined by the weights used in the cbcrc algorithm. 
In particular we explain how the weights in the cbcrc algorithm determine the shape 
of the search space of the generating vectors. In Section [4] we provide some numerical 
examples to illustrate that in certain instances the cbc2c algorithm is beneficial. 

2 The expectation value and standard deviation of 
the square worst-case error with random weights 

We assume that the weights j u , 7^ u C S, are nonnegative, independent random 
variables with a given mean and variance. Let E denote the expectation value and Var the 
variance. For any 7^ u C S there are numbers j u , a u > such that for all 7^ u, u' C S 
the following properties hold: 

• lu > 0, 

• E(7„) = %, 
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• E(7ti7«0 = E(7„)E(7 U /) for u ^ u', 

• Var(7„) = a 2 ; 

Note that we have E(t^) = 7« + o\- Let 7 = (jj^ucs and cr = (a u )^ u cs- 

We point out that the conditions for product weights need to be dealt with carefully. 
Assume that 7„ = rLeu^ f° r some nonnegative real numbers Zi. Further assume that 
E(zi) = Then E(7„) = Yl ieu Zi and hence 7 y u = Yl ieu Zi. Assume that Var(zi) = yf. 
Then for 0^u,t('C5 with v = u D u' we have 

E(7„7 U ') = E(JJ 2,; J J Zi) = ] [ Zi j ] Zj/ y*. 

Thus, using these assumption, we do not have E(7„7 U /) = E(7„)E(7 U /) in general, i.e., the 
weights are not independent (as is obvious from the definition 7 U = Yli^u z i)- The analysis 
for such weights is slightly different and is not considered here. Instead, if the weights are 
of product form 7„ = Yl ieu Zi, we assume that 7 y u = Yl ieu Zi and that a u = Yl ieu Vi for all 
7^ u C S for some numbers yi,Z{ > 0. 

If the weights 7„ are uniformly distributed in the interval [% + 5 U , 7 U — where < 
K < lui tri en the expectation value is 7 U > and E(7 2 ) = ^ f^ U -g™ it = l\ + 
thus the variance is 7^ = 5„/3. Note that 7^ — <5 U > 0, which ensures that the weights are 
always non-negative. In the following we do not assume that the weights are uniformly 
distributed, in fact, it is more interesting to assume a different distribution which also 
allows weights much larger than 7^ + 5 U . 

In the following the expectation E, the variance Var and the standard deviation Std 
are always taken with respect to the random variables 7 M . The expectation value is now 

1 "zl 

E(e 2 (U J ,P N (g)))= 7 u ^En^({^W) = e 2 (^,P,v^)). (2) 

Thus, current construction algorithms [El El [221 123] can be viewed as finding quadrature 
rules for which the expected value is small. Here we aim at finding quadrature rules for 
which, additionally, the variance is small. 

We point out that there is a difference between j u 'very small' and j u = 0. The 
restriction 7« > implies that if 7^ = then a u = 0. Thus if one constructs a lattice rule 
with error criterion e 2 (%y, P^(g)) where j u = for some 7^ u C S, then this means 
that the random variable 7 U = with probability 1. Thus setting 7 U = means that one 
knows that 7 U = and the associated ANOVA term f u = 0. Lattice rules constructed 
using such weights do not have any guarantee that positive weights 7« > will yield a 
good result. To illustrate, consider the two-dimensional example where 71 = 1, 72 = 1 
and 7i j2 = and the lattice rule has generating vector g = (1, 1). This lattice rule works 
well in this case but not if the weight 7^2 changes to 1, say; see also |20j. 

Using some elementary properties of the variance we obtain 

V^{e\U^P N {g)))= (^J^YlB 2 ({n gi /N})\ 
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and the standard deviation is given by 



Std(e 2 (H^,P N (g))) = JVaT(e 2 (n^P N (g))) 



\ 



£ f^EIWowiv})) 



Lattice rules for which the variance Var(e 2 ('H- ) ,, P^(g))) is small are less sensitive to 
changes of the weights j u . Notice that the variance is difficult to compute in general 
in high dimensions since it involves a sum over all subsets of {1, ... , s} for which o~ u > 
and which, in general, cannot easily be simplified to a formula which can be computed 
quickly, even in the case where the weights are of product form. 

Notice that ^ is the one-norm of the vector consisting of the error of the projections 
weighted by the expectation values of the weights, whereas the standard deviation is the 
two-norm of the vector consisting of the error of the projections weighted by the variance 
of the weights. 

Let fi be a probability measure on the weights (7 u )0^ M c,s- We now use the one-sided 
Chebyshev inequality which states that for a random variable X with probability measure 
Pr, expectation E(X) and standard deviation Std(X), we have for any c > that 

Pr(X -E(X) > cStd(X)) > 



1 + c 2 



Thus we obtain the following result. 
Lemma 1. For any c > we have 



fi ( 7 : e 2 (H^P N (g)) < e 2 {n^P N {g)) + c^d{e 2 {n^P N {g)))) > 

As noted above, the standard deviation is in general difficult to compute, however, 
using Jensen's inequality we have 



Std(e 2 (?C,iM<?))) 



N-l 



CS \ n=0 i£u 

N-l 



^ E °u^Y.Xl B ^ n 9i/ N })=z 2 {n (T ,P N {g)). 



Thus the square worst-case error with the variances as weights is an upper bound on the 
standard deviation. This upper bound can easily be computed (for instance for variances 
of product form). 

Assume now that if for some 7^ u C S we have a u = 0, then also 7^ = and therefore 
7 U = 0. Using Holder's inequality we have 

/ N-l 

e 2 (U~ f ,P N (g))= 7 ^EIl^^M) 



&u>0 



n=0 i£u 
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<Std(e 2 (W„,P K ( 9 ))) 



\ 0^«' 



7 2 
/ ii 



cs 
>o 



Thus a small standard deviation implies a small expected error. Combining the last two 
inequalities we obtain the following result. 

Lemma 2. Assume that if for some I ^ m C 5 we have a u = 0, then also 7« = and 
therefore 7 y u = 0. Then we have 



e 2 {H^P N {g)) < Std(e 2 CH CT , P N (g))) 



E ^<e 2 (H tT ,P N (g)) i 



E " 



^ o-l' 



If the weights cr are decaying such that the upper bound is independent of the di- 
mension, i.e. strong tractability (see [T^ [T5l 125] ) holds, then for any weights 7 such 

2 

that the expression XLcn m<oo ^ < 00 one a l so obtains a bound which is independent of 
the dimension, i.e. one has strong tractability. Further, if e 2 ('H (T , P/v(flO) satisfies strong 
tractability, then also the standard deviation is bounded. 

Combining the last two lemmas we obtain the following corollary. 

Corollary 1. Assume that if for some 7^ u C S we have a u = 0, then also 7„ = and 
therefore 7^ = 0. Then for any c > we have 



( 



( 



1 :e 2 {U^P N {g))<e 2 {U (T ,P N {g)) 



\ 



\ 



E 

\ <?u> 



,>0 



\ 



/ 



\ 



/ 



> 



1 + 6 



2 ' 



This result can be viewed as a robustness result with respect to weights. If one 
constructs a lattice rule using cr = (a u ) as weights, then for a set of weights taken from 
a set with measure at least jt^, the error is bounded by 



/ 



e 2 {U a ,P N {g)) 



\ 



E ^uM + C 



iC5 
4 >0 



There is one notable exception, namely, if <7 U = for some f) / u C 5, then no 
robustness with respect to the projection on u can be obtained (as was also illustrated 
above). 

In the following we show that one can do better by taking the robustness into account 
in the construction of the lattice rule itself. 



3 A construction of robust lattice rules 

In this section we generalize the component-by-component algorithm from [H [21]. 
repeat some facts from the fast cbc algorithm of Nuyens and Cools [161 E]- 



We 
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3.1 The fast Fourier transform method 



Nuyens and Cools (TBI E] have shown how to use the fast Fourier transform to reduce the 
computation time of the component-by-component algorithm. Because of the importance 
of these ideas we repeat them here (as is well understood, we see below that the algorithm 
of Nuyens and Cools actually calculates slightly more, which is important for the cbcrc 
algorithm below). 

For simplicity of exposition we assume product weights ^ u = Ylteu^i anc ^ that ^ 
is a prime number. For more general cases see [TBI Ej- Assume that the coordinates 
gl, . . . , g*_ x G {1, . . . , N — 1} are already fixed. Then we write the error criterion for 
j > 1 in the form 



e 2 (H- n P N ((g* 1 ,...,g*_ 1 ,z))) 

N-lj-l 

= - 1+ nT,U( 1+ Wing*/N})) (1 + %B 2 {{nz/N})) 



N 

n=0 i=l 



1 N-lj-l , N-lj-l 

■l + ^EII^ liMinQi/N})) + - II C 1 + %B 2 {{ngl/N})) B 2 {{nz/N}) 



n=0 i=l n=0 i=l 



1 ^ 

e\U^ P N ((g*, . . . , gU))) + - \{ (1 + %/&) /6 

i=l 



N-lj-l 



+ ^ E II C 1 + liMingl/N})) B 2 ({nz/N}). 



N 

n=l i=l 

Since the components gl,...,g*_ 1 are fixed, the value e 2 ("H-y, Pjv((<7i, • • • , Pj-i))) an d 
4 nti (1 + 7i/6) /6 does not depend on z and can therefore be ignored. Thus it suf- 
fices to calculate 

iv-i j-i 

#(z) = II ( X + liM{ngt/N})) B 2 ({nz/N}) for 1 < z < N - 1. 

n=l i=l 

We define the matrix 

O = (-B 2 ({712/iV}) ) !< nj;K <jv-i 

and the vector 

(j— 1 3-1 
n(i + 7^ 2 ({o^/iv})), , na + 7^ 2 ({(iv - ik*/iv})) 
i=l i=l 

Then 

(M/(l),...,vI/(iV-l)) T =fip. 

The matrix Q has some structure which allows one to use the fast Fourier transform. Let 
1 < v < N be a primitive element in the finite field Z^r = {0, 1, . . . , N — 1} of prime order 
N. Note that the multiplicative inverse v _1 is then also a primitive element. We define 
the permutation matrix IT(t>) = (irk i{ v ))i<k,i<N-i by 



1 if k = v l (mod N), 
otherwise. 
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Note that n(t>)IT(t>) T = J, the identity matrix. Let C = (ck,i)i<k,i<N-i be defined by 

C = U{v) T QU{v- 1 ), 

hence 

JV-l 

c k ,i = 7r Utk (v)B 2 ({uw/N}) 7 r W}l (v- 1 ) = B 2 {{v k ~ l / N}) . 

u,w=l 

The matrix C is therefore circulant. Let F/v-i = (N — i) -1 ^ 2 (fk,i)o<k,i<N-2 be the Fourier 
matrix of order N — 1 where fk,i = e 2nlkl ^ N . Then D = F^-xCF^^ is a diagonal matrix. 
Thus we have 

n = n^cn^-y = n(u)F J ^i 1J DF iV _ 1 n(?;- 1 ) T . 

Consider now the matrix- vector multiplication Qp. Multiplying a vector with the permu- 
tation matrices IT(t>), n(t> _1 ) T takes O(N) operations, the matrix vector-multiplication 
with the matrices F N _ ly F^\ can be carried out in O(NlogN) operations using the 
fast Fourier transform. Multiplying the diagonal matrix D with a vector takes O(N) 
operations. Thus the matrix-vector multiplication Qp can be carried out in 0(N log N) 
operations. For more details see [T8] . 

Notice that the fast matrix vector multiplication directly yields the whole vector 
. . . ,^(N — 1)) T . This vector can be ordered (using a sorting algorithm) to ob- 
tain ^{zi) < ^(z2) < • • • < ^(zn-i)- This can be done in O(NlogN) operations. Thus, 
by the above arguments, we can compute Zx, . . . , Zjv-i such that 

e 2 (H 7 , P N ((gl . . . , <?*_!, Zl ))) < ■ ■ ■< e 2 (H 7 , P N ((jg*, . . . , g*_ v z N ^))) 

in 0(N log N) operations. 



3.2 The fast cbc2c algorithm 

In the previous section we have shown some robustness of lattice rules which are con- 
structed for a given set of weights. In this section we modify the fast cbc algorithm jX6|, H7] 
to construct lattice rules for which, simultaneously, E(e 2 (H 7 , P/v (<?))) = e 2 (H 7 , P/v(</)) 
and Std(e 2 (% 7 , P/v(gr))) are small. Since the standard deviation Std(e 2 (H 7 , P/v(g))) is 
in general difficult to compute, we use e 2 (% cr , P/v(flO) as criterion instead. This has the 
additional advantage that the roles of 7 and a are interchangeable. 

Throughout the paper let k denote the number of distinct prime factors of the integer 
N > 2. We use jU Theorem 3], which states that for any c > 1, the proportion of 
generating vectors g E {1, . . . , N — 1} S which satisfy 

e'CH^PM) < E 7y r (2^- 2 C(2/r))H I for all 1 < r < 2, (3) 

where ((r) = YlT=i k~ r * s the Riemann zeta function, is bigger than 1 — c , i.e. there 
are more than (N — l) s (l — c -1 ) generating vectors g G {1, . . . , N — 1} S which satisfy the 
above bound (see also [19] for other criteria and bounds when N is not prime). Further, 
(U Theorem 10] states that a generating vector g* = (g\, . . . ,g*) e {1, . . . , iV — 1} which 
satisfies ^ can be found component-by-component. Thus we obtain the following result 
which follows from the fact that the intersection of two sets with measure bigger than 
1 — Ci 1 and 1 — C2 , where ci, C2 > 1 are such that 1 — c^ 1 + 1 — c 2 1 > 1, is non-empty. 
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Algorithm 1 (The fast component-by-component two criteria (fast cbc2c algorithm). 

Given: natural numbers N, s, nonnegative real numbers j u , o u for all ^ u C S; 1 < 
c i < 2, 2 < c 2 < 00 swc/i t/iat cj 1 + c 2 l = 1. 

• Setgl = 1; 

• For j = 2, . . . , s do the following: 

— Hard constraint 

Let K 1 = min{L(iV - 1)(1 - q 1 )] +1,N - 1}. Find the set of integers A = 
{zi, . . . , z Kl } C {1, . . . , iV — 1} which satisfies: 

e\H^ P N ((gt, gU^i))) < e2 (^7, PnM, <£-i, *))) 
/or all 1 < « < K\ and z E {1, . . . , N— 1}\A (using the fast Fourier transform) . 

— Soft constraint 

Let K 2 = min{[(iV - 1)(1 - c 2 l )\ + 1, N - 1}. Fznd i/ie set of integers B = 
{yi, . . . , 2/x 2 } Q {1j • • • ) ^ — 1} which satisfies: 

e\U a , P N ((gl, g*_ v &))) < e 2 ^, P JV ((^, . . . , y))) 

/or alll < i < K 2 andy G {1, . . . , iV— (using the fast Fourier transform). 

— Choose g* e An B which minimizes e 2 ((H^, P]y((gl, g*_ 1 ,w))) as a function 
ofw. 

• Return g* = (0*, . . . ,g* 8 ). 
Some comments are in order: 

(i) We have K x + K 2 > (N - 1)(1 - c\ 1 + 1 - c 2 l ) = N - 1. Thus the set AnB is not 
empty. Further note that the vector g* found by Algorithm [l] satisfies the bounds 
in Theorem [l] (see also jU Theorem 10]). 

(ii) The algorithm is basically symmetric in the constraints, but, by choosing 1 < c\ < 2, 
the first constraint is at least as hard to satisfy as the second one, since the upper 
bound is lower. 

(iii) We have biased the algorithm towards the hard constraint. Instead of choosing the 
value g* d G A fl B which minimizes e 2 ((H^, Pn({9*, ■ ■ ■ > 9d-n w ))) actually any value 
in the set An B could be chosen. The results still apply in this case. 

(iv) The classical cbc algorithm corresponds to the special case where c\ — 1 and c 2 = 00. 
Further, the classical cbc algorithm can also be obtained by choosing 7 = cr (in 
which case the choice of 1 < Cx, c 2 < 00 is irrelevant). 

(v) The fast cbc algorithm can be used to calculate the values e 2 ((H^, P/v((oi, • • • , 9*d-n z ))) 
and e 2 {H a , P N ((gl, . . . , g}_i,y))) for all y, z G {1, . . . , N - 1} very efficiently The 
values need to be sorted and then one needs to choose a value in A fl B. The 
main complexity is calculating the worst-case errors, hence the number of opera- 
tions needed for the algorithm is the same as that for the fast cbc construction 
[T6l [T7] (but with a larger constant since we have two worst-case errors). Thus one 
has a fast cbc2c algorithm. 
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From Theorem 10] we obtain the following result concerning the cbc2c algorithm. 

Theorem 1. Let N be an integer and 1 < ci, c-i < oo such that c\ x + 1 = 1. Then the 
generating vector g* 6 {1, . . . , N — 1} S constructed by the cbc2c algorithm satisfies 

e\H^ P N (g*)) < I | E 7y T (2 f %- 2 C(2/r))H I for all 1 < r < 2, and 

e 2 CH a , P N (g*))) < I | E ^ /T (2 K vr- 2 C(2/r))l«l I for alll<r<2. 
For c> and 1 < r, r' < 2 let 



T 



B, A c,r,r>N)=l% £ ^(2^- 2 C(2/r))N +c | £ ^'(2^(2//))'"' 

(4) 

We obtain the following corollary from Lemma [T] and Theorem [TJ 

Corollary 2. T/ie generating vector g* = (g^, . . . ,g*) constructed by the cbc2c algorithm 
satisfies 

c 2 

H (7 : e 2 CH 7 , P N (g*)) < B^(c, r, r' , N) for all 1 < r, r' < 2) > — — for any c > 0. 

By choosing 7 = cr, Algorithm [T] can be simplified to the classical cbc algorithm. Thus 
Corollary [2j with 7 = cr, applies to the classical fast cbc algorithm. However, in this case, 
if 7 U is small then also a u = 7 j u is small and thus the lattice rule may be sensitive to 
changes in the projection u. Unfortunately Corollary [2] does not give any information 
about the set of the weights which satisfy the condition. We study this topic in the 
following more general setting of the cbc algorithm with r constraints. 



3.3 The fast cbcrc algorithm 

In this subsection, instead of two constraints we study a cbc algorithm using 1 < r < 2 s — 1 
constraints. In this case one needs r sets of weights 7^, . . . ,7^ which are linearly in- 
dependent in ]R 2S_1 and 1 < ci, . . . , c r < 00 such that c^ 1 + • • • + c" 1 = 1. As we will 
see below, adding a vector of weights 7' which is a linear combination of the weights 
7W, . . . ,7^-* does not add a new constraint, since any vector satisfying the first r con- 
straints automatically satisfies the constraint using the weight 7'. 

Algorithm 2 (The fast component-by-component r criteria (fast cbcrc) algorithm). 

Given: natural numbers N,s and 1 < r < 2 s — 1, nonnegative real vectors 7^ = 
(7! )0^ u cs, ■ ■ ■ )7^ — (lu^)<d^ucs which are linearly independent in IR 2S_1 ; positive num- 
bers Ci, . . . , c r 6 R U {00} which satisfy c r > ■ ■ ■ > C\ > 1 and + ■ ■ ■ + c" 1 = 1. 

• Set K w = min{L(A- 1)(1 - c^ 1 )] + 1, N - 1} for 1 < w < r. 

• Set g\ = 1; 
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• For j = 2, . . . , s do the following: 

— For w = 1, . . . , r do the following: 

Find the set of integers A w = {z[ w \ . . . , z^} C {1, . . . ,N — 1} which satisfies: 

e 2 (U^ h P N ((gl g*^, z^))) < e 2 (^ 7 w, P N ((gl . . . , <£_!,*))) 

for all 1 < i < K w and z £ {1, . . . , iV — 1} \ A w using the fast algorithm 
described above. 

- Choose g* d £ V\x<w<r^w which minimizes e 2 (H 7 (i), P N ((g*, . . . , g* d _ x , z))) as a 
function of z. 

• Return g* = (gt,...,g*). 

The considerations above imply therefore that the construction cost of the fast cbcrc 
algorithm is 0(rsN(log N)) operations using 0(rN) storage (note that the intersection 
step can be done by sorting the elements in Ai,...,A r first, which takes 0(N log N) 
operations). 

From [4, Theorem 10] we also obtain a generalization of Theorem [T] which applies to 
the cbcrc algorithm. 

Theorem 2. Let N > 2 and r > 1 be integers and 1 < Ci,...,c r < oo such that 
c^ 1 + • • • + c" 1 = 1. Then the generating vector g* £ {1, . . . , N — 1} S constructed by the 
cbcrc algorithm using the weights . . . ,7^ satisfies 

e 2 (n^ h P N (g*)) < I ^ E (7^) 1/T (2 K 7r- 2 C(2/r))H I for all 1 < r < 2 
and all 1 < w < r . 



3.4 The geometry of the cbcrc algorithm 

For a £E £ IR 2S_1 we write x < y if Sj < yi for all 1 < i < 2 s , where x = (x±, . . . ,X2*-i) 
and y = (yi, . . . , ?/2 s -i)- Similarly we use the symbols <, >, >. 
For z = (zi, . . . , Z2°-i) £ IR 28 " 1 and e > we define the simplex 

T(z,e) = {yeR 2S - 1 :y>0,yz<e}. 

This simplex has vertices and (0, e/zi), 1 < % < 2 s , which stands for the vector whose z'th 
component is for i ^ i' and whose ith component is z$. If Zi = for some component 
i, then the simplex is degenerate and we consider the projection of the set onto those 
components which are nonzero (which is then a nondegenerate simplex). 
Let g = (gi, ...,g s ),g u = (gi) ieu and 

N-l 

e u = e u (g u ) = - ^ J] B 2 ({n 9l /N}) for ^ u C S, 

n=0 i€u 

e = e(g) = (e u ) 9 ^ u c S , 
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X 7 , c = X 7iC (iV) = M <2 ± £ (7«) 1/r (2 7 r- 2 C(2/r)) 

\ 0^«C5 

Notice that since e u is the worst-case error of integration in a reproducing kernel Hilbert 
space, we have e u > and therefore e > 0. 

Theorem [2] implies that the cbcrc algorithm now chooses the generating vector g e 
{1,...,JV- 1} S such that 

7 H . e := (7 H) . (Cu) ;= £ 7 M 6u < for ! < w < r 

Thus e lies in the intersection of the simplices r^w, X <„,) ) for 1 < w < r: 

eG f) r( 7 ^,X 7WiC J. (5) 

l<ui<r 

Geometrically this means that e lies in a convex r-polytope given by the intersection of r 
simplices. The weights 7^ change the shape of the simplices, whereas the values X 7 (») Cm 
change the size of the simplices. 

Compared to the cbc algorithm, the component-by-component r criteria (cbcrc) algo- 
rithm first increases the original simplex (by at most a factor of c r ) and then intersect it 
with other simplices, see Figure [T] This can be used to prevent e to be chosen too close 
to a vertex of the original simplex (this prevents e u from becoming too large for some 




Figure 1: The figure shows the search domains. The classical cbc algorithm searches in a 
simplex as indicated in the left picture. The cbc2c algorithm first increases the size of the 
simplex (second picture) and then intersects it with another simplex (third picture) to 
get the new search domain (fourth picture). This way extreme choices near the vertex of 
the original simplex (which corresponds to a large value of e u for some u) can be avoided. 



3.5 The geometry of the weights 

We now study the geometry of the weights for which the corresponding square worst-case 
error satisfies a certain bound. Since e is fixed once a generating vector g* is chosen, we 
consider now the set of weights 

r = {(7«W«c 5 : e^H^PNig*)) < ee 2 (H^Po)} , 
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where e > is a real number, P = and 



e 2 (n- f ,P )= inf 

|/||'H-v<l 



/(oj) da; 



[0,1] s 



is the initial error [25]. For our space we have e 2 (H 7 , Po) = !■ 

The square worst-case error e 2 (H 7 , P N (g*)) = e • 7 is a linear function of 7„. Thus T 
is a simplex in IR 2S_1 given by 

r = r( e ,e), 

which has vertices (0)^ u cs an d 

(O^e^e) for0^ M C5. 



Figure 2: The set of weights T. In the first instance, one value of e u (corresponding to the 
y-axis in the picture) is large, therefore the algorithm is sensitive to changes with respect 
to the projection u (the corresponding weight 7^ has to be small). A small change in 7^ 
could significantly increase the error. In the second case the algorithm is more robust 
since bigger changes in the second coordinate are allowed. 



We consider now the set of weights for which the cbcrc algorithm yields bounds. We 
have the following result. 

Theorem 3. Let g* be constructed by the cbcrc algorithm using the weights . . . , . 
Let 7 = Y^w=i ^7^' f or some Ai, . . . , A r > 0. Then it follows that 

r 

e 2 (H 1 ,P N (g*)) < y' j X w X^( W ) iCw . 

w=l 

Proof. By the cbcrc algorithm we have e 2 ("H 7 (™), P^(g*)) = e ■ 7^ < X y { W ) Cw for 1 < 
w < r. Thus we have 

r r r 

e\H^ P N (g*)) = e ■ 7 = £ X w e ■ 7^ = ^ \ w e 2 {U l{w) , P N (g*)) < ^ A^m^, 

W = l W = l UI = 1 

which shows the result. □ 

The above theorem shows that the cbcrc algorithm yields lattice rules which simulta- 
neously satisfy bounds for weights 7 G 1R 2S_1 taken from a subspace of IR 2 " -1 spanned by 

7W,..., 7 W el 2 '- 1 

The theorem above can be understood geometrically in the following way. Note that 
e G r(7 (u "), X 7 ( m ) jCm ) if and only if 7^ G r(e, X («) )Ch) ). We define the new weights 
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7 = eX ,L y^ w \ 1 < w < r for some real number e > 0. Then X^m = e 
for all 1 < w < r. Note that criterion in the cbcrc algorithm does not change by this 
normalization since the cbcrc algorithm yields exactly the same generating vector g* using 
the weights 7W, . . . , 7^ as it does for using the weights 7^, . . . , 7 . Thus 

e G r( 7 H x^J ^7 W G r(e,X 7W)C J ^ 7 W G r(e,e). 
Now ([5]) is equivalent to 

e G p| r( 7 ^, X 7WjCm ) ^ 7 (w) G r(e, e) for all 1 < w < r. 

l<w<r 

Therefore, the cbcrc algorithm ensures that g* is chosen such that 7 , . . . , 7^ all lie in 
the simplex T(e,e), i.e. 

7 (1) ,...,7 (r) Gr(e,e). 

In fact, the cbcrc algorithm finds, component-by-component, the smallest simplex T(e,e) 
which contains 7 , . . . , 7 . 

^ /|\ ('2 s — 1) 

If one chooses r = 2 s — 1 in the cbcrc algorithm, then the weights 7 V , . . . , 7 V ; G 
r(e,e) themselves are the vertices of a 2 s — 1-dimensional convex polytope which is con- 
tained in the simplex T(e,e) of the same dimension. However, if 1 < r < 2 s — 1, then 
the convex polytope spanned by 7 , . . . , 7^ is degenerate since it lies in a r-dimensional 
subspace. Thus using only r < 2 s — 1 weights 7*- 1 -*,. . . ,7^ does not fully control the 
shape of the simplex T(e,e). For the classical cbc algorithm only one vector of weights 
7 is used. The cbc construction then only ensures that 7 G T(e,e). In the numerical 
examples below we show that it is possible for the classical cbc construction to choose 
generating vectors which are not suitable for many other choices of weights. By adding 
additional constraints, the cbcrc algorithm can prevent such bad choices. 

The bound in Theorem [3] applies for all weights which lie in the linear subspace of M 2 " -1 
spanned by vectors 7^, . . . ,7^. In particular, if one uses the cbc(2 s — l)c algorithm, 
then one can obtain a bound for any choice of weights. However, in higher dimensions s 
this is currently problematic for two reasons: the computational cost is exponential in the 
dimension; the second problem is that one would have to choose Ci, . . . , C2«_i such that 
c^ 1 H — • + o%t_x = 1- For instance, the choice Cj = 2 s — 1, 1 < j < s, yields a factor in the 
upper bound which grows exponentially with the dimension. For lower dimensions this is 
feasible though and hence can be useful in applications with low truncation dimension. 

Consider now 7 = J2w=i > where Ai, . . . , A r > and Ai + • • • + A r < 1, i.e., 7 

lies in the convex polytope with vertices 0, 7^, . . . , 7^. Then, by Theorem [3J we have 

e 2 CH^P N (g*)) <£■ 
We summarize the results in the following corollary. 

Corollary 3. Let 1 < r < 2 s — 1. Let c\, . . . , c r > 1 be given such that c^ 1 + • • ■ + c~ 1 = 1. 

Let e > be a real number. Let 7 , . . . , 7^ be weights which are normalized such 
that X~( W ) Cw = £ for all 1 < w < r. Let 7 = J2w=i ^wT' 1 "' where Ai, . . . , A r > and 
Ai + • • • + A r < 1. Let g* be constructed by the cbcrc algorithm based on the weights 
7 (1) ,...,7 {r) . Then 

7 (1) ,...,7 (r) Gr(e,e) 

and 

e 2 (H 7 ,P^(y*)) <e. 
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4 Numerical results 



To illustrate the ideas in the paper we chose some instructive examples. We tested the 
algorithm with r = 2. The computation time for the cbc2c algorithm is between 2 and 
3 times the computation time of the cbclc algorithm. This was observed for a variety of 
choices for c\ and C2 and values n. 

Further we tested the component-by- component algorithm with fast decaying weights. 
We used s = 100, product weights 7^ = njeuT? w hh tfp = 1, 7^ = lO - ^ and 

7^ = Further we chose c\ = C2 = 2. The cbclc algorithm using the weights 
returns the same components g* = g*, for j,j' > 15. This choice of generating vector 
would be bad for slow decaying weights. On the other hand, the cbc2c does not return 
any repeated components, which is prevented by the second constraint. The results are 
presented in Table [T] It shows that the cbc2c algorithm yields approximately the same 
results as the cbclc algorithm constructed for the right weight, but can do significantly 
better if the lattice rule constructed by the cbclc algorithm is used for different weights, 
as can be seen in Table [T] The results in Table [2] are for a different choice of weights and 
are similar. 

The choice of weights 7^ and jj 2 ^ in Table [l] and Table |i] are quite different from 
each other. Numerical tests for examples where the weights are more similar than in 
the examples shown, for instance 7- = j _1 and t ? - 2 ' ) = j~ 4 (or even j^. = 2~ J ), yield 
numerical results which are quite similar, indicating that there are not many (bad) outliers. 
Table [3] and [4] show numerical results with randomly chosen weights, again showing that 
the lattice rules constructed by the cbc and cbc2c algorithms perform well except in the 
case where the cbc construction is based on fast decaying weights. 
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N 


251 


509 


1019 


2039 


4079 


cbc2c 


e( 7 (1) ) 


1.4044e+02 


9.8623e+01 


6.9702e+01 


4.9275e+01 


3.4838e+01 


cbc( 7 «) 


e( 7 (1) ) 


1.4044e+02 


9.8623e+01 


6.9702e+01 


4.9274e+01 


3.4838e+01 


cbc( 7 ( 2 )) 


e(7 (1) ) 


2.4075e+02 


2.1790e+02 


1.9762e+02 


1.9481e+02 


1.8137e+02 


cbc2c 


e(7 (2) ) 


5.4897e-04 


2.7128e-04 


1.3568e-04 


6.7927e-05 


3.3965e-05 


cbc( 7 «) 


e( 7 (2) ) 


5.4897e-04 


2.7128e-04 


1.3568e-04 


6.7930e-05 


3.3966e-05 


cbc( 7 ( 2 )) 


e(7 (2) ) 


5.4882e-04 


2.7113e-04 


1.3558e-04 


6.7892e-05 


3.3954e-05 


cbc2c 


e ( 7 (3)) 


3.197le-02 


1.9872e-02 


1.2057e-02 


7.7449e-03 


4.9349e-03 


cbc( 7 ( 1 ') 


e ( 7 (3)) 


3.1653e-02 


2.0250e-02 


1.3100e-02 


8.5062e-03 


4.9972e-03 


cbc( 7 ( 2 )) 


e ( 7 (3)) 


l.5597e-01 


1.4616e-01 


1.3816e-01 


1.3813e-01 


1.3130e-01 


N 


8161 


16319 


32633 


65267 


130531 


cbc2c 


e( 7 W) 


r\ a r* c\r\ i r\ i 

2.4629e+01 


1.7417e+01 


1.2316e+01 


8.7088e+00 


6.1579e+00 


cbc(V 1 )) 


e( 7 (1) ) 


2.4629e+01 


1.7417e+01 


1.2316e+01 


8.7087e+00 


6.1579e+00 


Cbc( 7 (2)) 


e( 7 (1) ) 


1.6777e+02 


1.6737e+02 


1.5575e+02 


1.5567e+02 


1.5563e+02 


cbc2c 


e( 7 (2) ) 


1.7023e-05 


8.5236e-06 


4.2695e-06 


2.1370e-06 


1.0753e-06 


cbc( 7 «) 


e( 7 (2) ) 


1.7023e-05 


8.5236e-06 


4.2690e-06 


2.1370e-06 


1.0721e-06 


cbc( 7 ( 2 )) 


e( 7 (2) ) 


1.7006e-05 


8.5111e-06 


4.2631e-06 


2.1351e-06 


1.0683e-06 


cbc2c: e( 7 (3) ) 


3.0911e-03 


2.0308e-03 


1.2551e-03 


7.9994e-04 


5.2220e-04 


cbc(V 1 ') 


e ( 7 (3)) 


3.2965e-03 


2.1254e-03 


1.3246e-03 


8.5575e-04 


5.8758e-04 


cbc( 7 ( 2 )) 


e ( 7 (3)) 


1.2522e-01 


1.2521e-01 


1.1979e-01 


1.1979e-01 


1.1979e-01 


N 


261061 


522127 


1044257 


2088511 


4177051 


cbc2c 


e( 7 (1) ) 


4.3542e+00 


3.0787e+00 


2.1769e+00 


1.5392e+00 


1.0883e+00 


cbc(V 1 )) 


e( 7 (1) ) 


4.3542e+00 


3.0787e+00 


2.1769e+00 


1.5392e+00 


1.0883e+00 


cbc( 7 ( 2 )) 


e( 7 (1) ) 


1.3495e+02 


1.1735e+02 


1.2566e+02 


1.1702e+02 


1.1701e+02 


cbc2c 


e( 7 (2) ) 


5.3706e-07 


2.7758e-07 


1.3899e-07 


6.3220e-08 


0.0000e+00 


cbc(7«) 


e( 7 (2) ) 


5.3686e-07 


2.5938e-07 


1.3899e-07 


0.0000e+00 


7.3000e-08 


cbc( 7 ( 2 )) 


e( 7 (2) ) 


5.3541e-07 


2.5981e-07 


1.3411e-07 


1.1151e-07 


1.6255e-07 


cbc2c 


e ( 7 (3)) 


3.2756e-04 


2.1752e-04 


1.4107e-04 


8.6973e-05 


5.7966e-05 


cbc(7«) 


e ( 7 (3)) 


4.2678e-04 


2.1369e-04 


1.3647e-04 


9.4174e-05 


5.7208e-05 


cbc( 7 ( 2 )) 


e ( 7 (3)) 


1.1070e-01 


1.0250e-01 


1.0620e-01 


1.0234e-01 


1.0234e-01 



Table 1: This table shows the square worst-case errors using the cbc2c construction 
based on the weights 7^ and 7^, the cbc construction based on the weights 7^ and 
the cbc construction based on the weights 7^. Here, e( 7 ^) stands for the worst-case 
error e(H^(w), P?v(gr*)). We choose s = 100, product weights with 7^ = 1, = 10~ j , 
if =J-\ c 1 = c 2 = 2; 
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N 


251 


509 


1019 


2039 


4079 


cbc2c: e( 7 ilJ ) 
cbc( 7 «): e( 7 «) 
cbc( 7 ( 2 )): e^ 1 )) 


3.9170e-03 
3.9148e-03 
9.5094e-03 


2.0603e-03 
2.0587e-03 
7.6012e-03 


1.1145e-03 
1.1142e-03 
6.4146e-03 


6.0564e-04 
6.0547e-04 
6.3616e-03 


3.2079e-04 
3.2073e-04 
5.5514e-03 


cbczc: e(7 l ; ) 
cbc( 7 ( 1 )): e( 7 ( 2 )) 
cbc( 7 ( 2 )): e( 7 ( 2 )) 


5.4893e-04 
5.4893e-04 
5.4882e-04 


2.7113e-04 
2.7113e-04 
2.7113e-04 


1.3558e-04 
1.3558e-04 
1.3558e-04 


6.7914e-05 
6.7914e-05 
6.7892e-05 


3.3957e-05 
3.3957e-05 
3.3954e-05 


cbc2c: e( 7 w ) 
cbc( 7 W): e( 7 ( 3 )) 
cbc( 7 ( 2 )): e( 7 ( 3 )) 


5.3522e-03 
5.9244e-03 
3.0073e-02 


3.2454e-03 
2.9483e-03 
2.9957e-02 


2.1131e-03 
1.9003e-03 
2.9922e-02 


r\ An 1 1 f\ A 

9.98lle-04 
l.02l5e-03 
2.99l4e-02 


5.6878e-04 
5.9843e-04 
2.9906e-02 


N 


8161 


16319 


32633 


65267 


130531 


cbc2c: e( 7 ii; ) 
cbc^ 1 )): e( 7 «) 
cbc( 7 ( 2 )): e( 7 «) 


1.7664e-04 
1.7664e-04 
4.9127e-03 


9.5576e-05 
9.5576e-05 
4.9093e-03 


5.2562e-05 
5.2559e-05 
4.3999e-03 


2.86l6e-05 
2.86l6e-05 
4.3994e-03 


1.5674e-05 
1.5674e-05 
4.3993e-03 


cbc2c: e( 7 ^ 2 )) 
cbc( 7 W): e( 7 ( 2 )) 
cbc( 7 ( 2 )): e( 7 ( 2 )) 


1.7006e-05 
1.7006e-05 
1.7006e-05 


8.5147e-06 
8.5147e-06 
8.5111e-06 


4.2646e-06 
4.2646e-06 
4.2631e-06 


2.l34le-06 
2.l34le-06 
2.l35le-06 


1.0650e-06 
1.0650e-06 
1.0683e-06 


cbc2c: e( 7 (3) ) 
cbc( 7 «): e( 7 ( 3 )) 
cbc( 7 ( 2 )): e( 7 ( 3 )) 


3.5384e-04 
3.5384e-04 
2.9898e-02 


1.9131e-04 
1.9131e-04 
2.9898e-02 


1.1179e-04 
1.0812e-04 
2.9891e-02 


6.73l0e-05 
6.73l0e-05 
2.989le-02 


a a no 1 r\ v 

4.883le-05 
4.2383e-05 
2.989le-02 


N 


261061 


522127 


1044257 


208851 1 


4177051 


cbc2c 
cbc( 7 W) 
cbc( 7 ( 2 )) 


e( 7 (1) ) 
e( 7 (1) ) 
e(7 (1) ) 


8.6033e-06 
8.6019e-06 
3.6719e-03 


4.7039e-06 
4.7038e-06 
3.0833e-03 


2.6006e-06 
2.5996e-06 
3.3182e-03 


1.4241e-06 
1.4242e-06 
3.0560e-03 


8.02l8e-07 
8.0024e-07 
3.0560e-03 


cbc2c 
cbc( 7 «) 
cbc( 7 ( 2 )) 


e(7 (2) ) 
e(7 (2) ) 
e( 7 (2) ) 


5.3458e-07 
5.3458e-07 
5.3541e-07 


2.7028e-07 
2.7028e-07 
2.5981e-07 


1.1921e-07 
1.1921e-07 
1.3411e-07 


0.0000e+00 
0.0000e+00 
1.1151e-07 


0.0000e+00 
0.0000e+00 
l.6255e-07 


cbc2c 
cbc( 7 «) 
cbc( 7 ( 2 )) 


e ( 7 (3)) 
e ( 7 (3)) 
e ( 7 (3)) 


2.0703e-05 
2.3378e-05 
2.9876e-02 


1.2833e-05 
1.2218e-05 
2.9860e-02 


7.2400e-06 
7.6749e-06 
2.9868e-02 


5.4021e-06 
3.9504e-06 
2.9860e-02 


2.73l4e-06 
2.0520e-06 
2.9860e-02 



Table 2: This table shows the square worst-case errors using the cbc2c construction 
based on the weights 7^ and 7^, the cbc construction based on the weights 7^ and 
the cbc construction based on the weights 7^. Here, e( 7 ^) stands for the worst-case 
error e^K^.P^ig*)). We choose s = 100, product weights with 7^ = l, = j~ 2 , 
7f = (s-jr 2 ,c l = c 2 = 2; 
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N 


251 


509 


1019 


2039 


4079 


cbc2c: e(7^) 
cbc(7«): e(7«) 
cbc( 7 ( 2 )): e(7«) 


3.0814e-02 
2.9982e-02 
8.5356e-02 


1.8765e-02 
1.8448e-02 
7.9743e-02 


1.1539e-02 
1.1482e-02 
7.7126e-02 


7.0582e-03 
7.0420e-03 
6.7998e-02 


4.4014e-03 
4.3849e-03 
6.7529e-02 


cbc2c: e(7 (2) ) 
cbc( 7 ( 1 )): e( 7 ( 2 )) 
cbc( 7 ( 2 )): e( 7 ( 2 )) 


2.4617e-03 
2.4617e-03 
2.4416e-03 


1.2385e-03 
1.2385e-03 
1.2423e-03 


6.6254e-04 
6.6254e-04 
6.5820e-04 


3.5517e-04 
3.5517e-04 
3.4793e-04 


1.8239e-04 
1.8239e-04 
1.7957e-04 


cbc2c: e(7 w ) 
cbc( 7 ( 1 )): e( 7 ( 3 )) 
cbc( 7 ( 2 )): e( 7 ( 3 )) 


2.5477e+00 
2.5413e+00 
4.0922e+00 


1.7828e+00 
1.7826e+00 
4.7043e+00 


1.2581e+00 
1.2576e+00 
3.7287e+00 


8.8750e-01 
8.8690e-01 
3.1589e+00 


6.2621e-01 
6.2570e-01 
3.2877e+00 


N 


8161 


16319 


32633 


65267 


130531 


cbc2c: e(7 l J ) 
cbc( 7 ( 1 )): e(7«) 
cbc( 7 ( 2 )): e(7«) 


2.7229e-03 
2.7208e-03 
6.1343e-02 


1.6958e-03 
1.6957e-03 
5.5530e-02 


1.0601e-03 
1.0587e-03 
5.3283e-02 


6.6402e-04 
6.6282e-04 
5.2673e-02 


4.1363e-04 
4.1368e-04 
5.0934e-02 


cbc2c: e(7 l - 1 ) 
cbc( 7 ( 1 )): e( 7 ( 2 )) 
cbc( 7 ( 2 )): e( 7 ( 2 )) 


9.5776e-05 
9.5776e-05 
9.4743e-05 


4.9506e-05 
4.9506e-05 
4.9263e-05 


2.6222e-05 
2.6222e-05 
2.5759e-05 


1.3943e-05 
1.3943e-05 
1.3567e-05 


7.3677e-06 
7.3677e-06 
7.2127e-06 


i ~ / /q\\ 

cbc2c: e(7 (3) ) 
cbc( 7 ( 1 )): e( 7 ( 3 )) 
cbc( 7 ( 2 )): e( 7 ( 3 )) 


4.4091e-01 
4.4099e-01 
3.1957e+00 


3.1093e-01 
3.1075e-01 
2.6282e+00 


2.1929e-01 
2.1886e-01 
2.7940e+00 


1.5442e-01 
1.5431e-01 
2.5026e+00 


1.0883e-01 
1.0875e-01 
2.3333e+00 


N 


261061 


522127 


1044257 


2088511 


4177051 


cbc2c: e(jW) 
cbc(7«): e(7W) 
cbc( 7 ( 2 )): e(7«) 


2.5887e-04 
2.5874e-04 
4.6470e-02 


1.6206e-04 
1.6160e-04 
4.7789e-02 


1.0112e-04 
1.0101e-04 
4.6729e-02 


6.3253e-05 
6.3241e-05 
4.1838e-02 


3.9582e-05 
3.9495e-05 
4.5181e-02 


cbc2c: e(j^) 
cbc( 7 ( 1 )): e( 7 ( 2 )) 
cbc( 7 ( 2 )): e( 7 ( 2 )) 


3.8690e-06 
3.8690e-06 
3.7567e-06 


2.0481e-06 
2.0481e-06 
1.9802e-06 


1.0884e-06 
1.0884e-06 
1.0503e-06 


5.8400e-07 
5.8419e-07 
5.3995e-07 


2.8115e-07 
2.6781e-07 
3.0465e-07 


cbc2c: e(7 (3) ) 
cbc( 7 ( 1 )): e( 7 ( 3 )) 
cbc( 7 ( 2 )): e( 7 ( 3 )) 


7.6712e-02 
7.6574e-02 
2.2229e+00 


5.3901e-02 
5.4043e-02 
2.1986e+00 


3.7986e-02 
3.7979e-02 
2.3036e+00 


2.6700e-02 
2.6664e-02 
2.0369e+00 


1.8762e-02 
1.8764e-02 
2.1620e+00 



Table 3: This table shows the square worst-case errors using the cbc2c construction based 
on the weights 7^ and 7*- 2 - 1 , the cbc construction based on the weights 7*- 1 ) and the 
cbc construction based on the weights 7^. Here, e(j^) stands for the worst-case error 
e(H^(w),PN(g*))- We choose s = 100, product weights with 7^ = j" 1 , 7^ = 2~ j , 7^ is 
chosen randomly in [0, 1], c x = c 2 = 2; 
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N 


or 1 

251 


r nn 

509 


1 n 1 n 

1019 


ono n 

2039 


a n Tn 

4079 


cbc2c: 


e(7 (1) ) 


O A A 1 P.^ no 

z.441oe-Uo 


1 o/ioQn no 
l.z4zoe-Uo 


p. coon^ dA 
0.0ozUe-U4 


Q ^7fl7^ CiA 


1 7ncc7^ n/i 

1. / yo 1 e-U4 


cbc( 7 ( 1 )): 


e( 7 (1) ) 


O A A "1 r* f\ o 

2.4416e-03 


1.2423e-03 


t~> p* o o n n a 

6.5820e-04 


o A Tn o n A 

3.4793e-04 


1 Tn r* T n /I 

1.7957e-04 


cbc( 7 ( 2 )): 


e( 7 (1) ) 


O A O OO C\ O 

2.4832e-03 


1 o Tn r no 

1.2705e-03 


r* r* a nc\ n a 

6.6492e-04 


o r 1 1 no n 1 

3.6102e-04 


1 O TO i n A 

1.8731e-04 


cbc2c: 


e(7 (2) ) 


o /inno~ i nn 
o.4(Jyoe+UU 


o oni i nn 

z.oyl * e+UU 


1.6oooe+UU 


1 1 noo^ i nn 


O A Ol ni 


cbc( 7 ( 1 )): 


e( 7 (2) ) 


r or /i r i nn 

5.3545e+00 


/i i nrr i nn 

6.l955e+00 


a o o o "i i nn 

4.8321e+00 


o nnin i nn 

3.9970e+00 


/i o r* o r* i nn 

4.2535e+00 


cbc( 7 ( 2 )): 


e(7 (2) ) 


o /inro i nn 

3.4(J52e+(J(J 


o ooto i nn 

2.3878e+00 


1 /in i nn 

1.6849e+00 


i i oon i nn 

1.1889e+00 


o o o t /i n 1 

8.3874e-01 


cbc2c: 


e ( 7 (3)) 


1 A CiA A rt l no 

1.4U44e+Uz 


n ocoo^ i m 


o nvno~ i ni 

6.y / Uoe+Ul 


4.yz t oe+Ul 


o /i ooo„ i ni 

o.4oooe+Ul 


cbc^ 1 )): 


e ( 7 (3)) 


i A 'in i i no 

1.4761e+02 


1 on r o i no 

1.2953e+02 


o /i r o a i n 1 

8.4534e+01 


6.2yi6e+01 


r rnno i n 1 

5.5008e+01 


cbc( 7 ( 2 )): 


e ( 7 (3)) 


1 A C\ A A i no 

1.4044e+02 


n o po o i n 1 

9.8623e+01 


p nTnn i ni 

6.9702e+01 


/i noTr* i n 1 

4.9275e+01 


o a o o o i n 1 

3.4838e+01 


N 


8161 


16319 


OO (~* o o 

32633 


/" r" op t 

65267 


1 on r* o 1 

130531 


cbc2c: 


e(7 (1) ) 


y.4/4oe-Uo 


4.yzooe-Uo 


z.o /oye-Uo 


1 Qceo„ np; 


7 oi io„ nt? 


cbc^ 1 )): 


e(7 (1) ) 


n 1 t 1 o n r - 

9.4743e-05 


/i nor* o n r 

4.9263e-05 


o rTrn nr 

2.5759e-05 


1 o r* p T nr 

1.3567e-05 


t o 1 ot nn 

7.2127e-06 


cbc( 7 ( 2 )): 


e(7 (1) ) 


9.6866e-05 


r* 1 Tno nr 

5.1703e-05 


o p o a o n r* 

2.6248e-05 


1 /I r* n T n r* 

1.4507e-05 


T T a o n n z" 1 

7.7439e-06 


cbc2c: 


e(7 (2) ) 


aA n 1 

o.y4obe-Ul 


a 1 nno^ m 

4.iyyze-ui 


z.yo4 / e-Ui 


o nnoi ^ m 

z.uyzle-Ul 


1.4 fooe-Ul 


cbc( 7 «): 


e(7 (2) ) 


o r* a r* o inn 

3.6458e+00 


o o tp a i nn 

3.2764e+00 


o 1 nn i i nn 

3.1991e+00 


o o 1 n /i inn 

3.2104e+00 


o o i i /i i nn 

3.3114e+00 


cbc( 7 ( 2 )): 


e(7 (2) ) 


r ai nc\ n i 

5.9162e-01 


/i 1 to i n 1 

4.1731e-01 


o n 1 o i n 1 

2.9421e-01 


o nToo ni 

2.0732e-01 


1 /i r*n r n 1 

1.4605e-01 


cbc2c: 


e ( 7 (3)) 


z.4oo(Je+Ul 


1 T A 1 i n 1 


1 oo i v~ i ni 

l.zol /e+Ul 


q 7nno^ i nn 


el i r-oo„ i nn 


cbc( 7 «): 


e ( 7 (3)) 


a a nro i A 1 

4.1(J63e+(Jl 


o 1 on t i n 1 

3.12y7e+01 


o to r n i n 1 

2.7350e+01 


o roon i ni 

2.522ye+01 


o oor ,i i n 1 

2.2864e+01 


cbc( 7 ( 2 )): 


e ( 7 (3)) 


o /i nnn i n 1 

2.4629e+01 


1 i a a i i n 1 

1.7417e+01 


1 o ' ) i < ■ i ni 

1.2316e+01 


o Tnoo i nn 

8.7088e+00 


p iron i nn 

6.1580e+00 


N 


zolUol 


C001 <)7 


1 CiA 

1U44Z0 ( 




a 1 77n 1 
41 f /Ui3l 


cbc2c: 


e(7 (1) ) 


3.7566e-06 


1.9798e-06 


1.0551e-06 


5.4057e-07 


2.9352e-07 


cbc( 7 «): 


e(7 (1) ) 


3.7567e-06 


1.9802e-06 


1.0503e-06 


5.3995e-07 


3.0465e-07 


cbc( 7 ( 2 )): 


e(7 (1) ) 


4.3237e-06 


2.1728e-06 


1.2567e-06 


6.4075e-07 


3.1434e-07 


cbc2c: 


e(7 (2) ) 


1.0412e-01 


7.3434e-02 


5.1821e-02 


3.6553e-02 


2.5707e-02 


cbc( 7 «): 


e(7 (2) ) 


3.0074e+00 


2.9431e+00 


3.1018e+00 


2.7223e+00 


2.8919e+00 


cbc( 7 ( 2 )): 


e(7 (2) ) 


1.0284e-01 


7.2401e-02 


5.0932e-02 


3.5823e-02 


2.5188e-02 


cbc2c: 


e ( 7 (3)) 


4.3546e+00 


3.0791e+00 


2.1772e+00 


1.5395e+00 


1.0886e+00 


cbc^ 1 )): 


e ( 7 (3)) 


1.8153e+01 


1.9073e+01 


1.8373e+01 


1.4183e+01 


1.6737e+01 


cbc( 7 ( 2 )): 


e ( 7 (3)) 


4.3543e+00 


3.0788e+00 


2.1770e+00 


1.5393e+00 


1.0884e+00 



Table 4: This table shows the square worst-case errors using the cbc2c construction based 
on the weights and the cbc construction based on the weights and the 
cbc construction based on the weights 7^. Here, e(7^) stands for the worst-case error 
e(H^(w),PN{g*))- We choose s = 100, product weights with 7]^ = 2~ j , rfp is chosen 
randomly, 7j 3) = 1, c x — c 2 = 2; 
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